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Abstract 

In  order  to  describe  the  dynamic  behaviour  of  a  fuel  cell  system  its  components  are  modelled  by  the  help  of  ID  dynamic  models  which 
are  implemented  in  Matlab  Simulink.  The  fuel  cell  system  consists  of  an  autothermal  gasoline  reformer  (ATR)  which  for  the  realisation  of 
a  high  system  efficiency  is  thermally  coupled  to  the  other  system  components  (gas  purification,  heat  exchangers).  Dynamic  simulations  of 
load  changes  show  that  the  dynamic  behaviour  of  such  a  system  is  primarily  dominated  by  the  response  times  of  the  liquid  water  flowing 
through  the  heat  exchangers,  the  volume  of  which  should  consequently  be  reduced  to  a  minimum  in  order  to  achieve  shorter  response 
times.  In  contrast,  the  dynamic  behaviour  of  the  reactors  is  not  critical.  The  composition  of  the  product  gas  at  the  ATR-outlet  is  however 
influenced  by  the  moisture  content  of  the  gas  at  the  ATR  inlet,  which  in  turn  is  negatively  influenced  for  a  short  transition  period  by  the 
residence  times  of  the  water  flowing  through  the  heat  exchangers  during  a  load  change. 

According  to  the  results  obtained,  the  system  is  able  to  adjust  to  load  changes  within  20  s  for  a  load  increase  (10-90%  of  full  load)  and 
within  3  s  for  a  load  decrease  (90-10%  of  full  load). 
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1.  Introduction 

Fuel  cell  systems  are  considered  to  have  a  high  potential  to 
increase  the  efficiency  of  modern  cars  [1-3].  In  order  to  make 
use  of  the  existing  infrastructure  and  to  use  the  high-energy 
density  of  gasoline  there’s  a  lot  of  research  going  on  in 
gasoline  reforming  for  on-board  hydrogen  production  in  fuel 
cell  vehicles  [4-6]. 

To  be  an  alternative  to  modern  combustion  engines,  fuel 
cell  systems  with  hydrogen  production  by  gasoline  reform¬ 
ing  have  to  meet  the  following  requirements: 

•  Fuel  cell  systems  have  to  enable  higher  efficiencies  than 
internal  combustion  engines  in  order  to  lower  the  fuel 
consumption  of  modern  cars. 

•  The  start-up  time  and  the  dynamic  behaviour  of  such  a 
vehicle  has  to  be  comparable  to  the  one  of  modern  cars. 

•  The  emissions  of  a  fuel  cell  vehicle  have  to  meet  the 
future  emission  regulations  regarding  CO,  HC,  NOx  and 
particles  (e.g.  SULEV,  Euro  V). 
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•  They  have  to  be  competitive  with  regard  to  volume,  weight 
and  costs. 

These  requirements  for  fuel  cell  systems  with  on-board 
gasoline  reformers  are  affected  in  particular  by  the  operat¬ 
ing  parameters  and  the  characteristics  of  the  reformer  whose 
efficiency  is  given  by  the  heating  value  of  the  produced 
hydrogen  in  relation  to  the  heating  value  of  the  fed  gasoline: 

_  m2  •  a//h2 

^/Reformer  —  T  .  u  • 

^Gasoline  *  A /7  Gasoline 

For  the  hydrogen  production  by  gasoline  reforming  most 
publishers  [7,8,4]  prefer  autothermal  reforming  concepts  to 
steam  reforming  and  partial  oxidation  because  they  enable 

•  a  high  hydrogen  yield  because  of  the  addition  of  water  to 
the  feed; 

•  minimisation  of  NOx  —  and  soot-production  by  the 
addition  of  water  and  the  low  reaction  temperatures 
(800-1000  °C); 

•  dynamic  operation  through  in-situ-provision  of  the  re¬ 
quired  energy  because  of  exothermal  reactions. 

The  operating  parameters  of  autothermal  reformers  are  the 
steam-to-carbon  ratio  S/C  which  brings  the  fed  water  into 
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Nomenclature 

CLy/ 

specific  catalyst  surface  area  (m2/m3) 

A 

cross-section  area  of  one  channel  (m2) 

cp,  g 

heat  capacity  gas  (J/mole  K) 

Cs 

heat  capacity  solid  (J/kg  K) 

A,  g 

diffusions  coefficient  of  component  i 
in  the  gas  phase  (mole/m  s) 

h 

enthalpy  of  water  (J/mole) 

h! 

enthalpy  of  liquid  water  at  boiling 
point  (J/mole) 

hc 

carbon  mole  flow  (mole/s) 

ho 

gas  mole  flow  (mole/s) 

^Gasoline 

gasoline  mole  flow  (mole/s) 

«h2 

hydrogen  mole  flow  (mole/s) 

^h2o 

water/steam  mole  flow  (mole/s) 

«o2 

oxygen  mole  flow  (mole/s) 

^02,stoich. 

oxygen  mole  flow  for  stoichiometric 
combustion  (mole/s) 

Rj 

reaction  rate  of  reaction  j  (mole/m2  s) 

S/C 

steam-to-carbon  ratio 

SR 

stoichiometric  ratio 

7g 

gas  temperature  (K) 

Ts 

solid  temperature  (K) 

u 

velocity  (m/s) 

gas  flow  rate  (m2/s) 

yi,  g 

mole  fraction  of  component  i  in  the 
gas  phase 

yi,s 

mole  fraction  of  component  i  at  the 
catalyst  surface  area 

as,g 

heat  transfer  coefficient  gas  — >►  solid 
(W/m2  K) 

Pi 

mass  transfer  coefficient  of 
component  i  (m/s) 

A/Zb 

enthalpy  for  vaporization  (J/mole) 

AF/Gasoline 

lower  heating  value  of  gasoline  (J/mole) 

A% 

lower  heating  value  of  hydrogen  (J/mole) 

A  Hj 

enthalpy  of  reaction  j  (J/mole) 

£ 

porosity  of  catalyst  structure 

^Reformer 

reformer  efficiency 

kg 

heat  conductivity  gas  (W/m  K) 

heat  conductivity  solid  (W/m  K) 

viJ 

stoichiometry  coefficient  of  component  i, 
reaction  j 

Pg 

gas  density  (mole/m2) 

Ps 

solid  density  (kg/m2) 

relationship  with  the  carbon  in  the  feed 

s/c  =  ^, 

he 

and  the  stoichiometric  ratio  SR  which  is  defined  as  the 
amount  of  oxygen  in  the  feed  divided  by  the  amount  of 


oxygen  necessary  for  complete  combustion: 


^02,stoich. 

In  addition  to  these  operating  parameters,  the  reactor  tem¬ 
perature  of  an  autothermal  reformer  is  influenced  by  the 
feed  temperature. 

In  order  to  achieve  high-system  efficiencies,  a  suitable 
system  configuration  was  developed  with  the  help  of  the 
Pinch  Point  method  and  steady- state  simulations  (Aspen 
Plus),  indicating  a  probable  system  efficiency  of  30%.  The 
fuel  cell  system,  whose  basic  design  is  shown  in  Fig.  1 
comprises  the  components  ATR  for  the  conversion  of  gaso¬ 
line  into  hydrogen,  high  temperature  shift  reactor  (HTS), 
low  temperature  shift  reactor  (LTS),  preferential  oxidation 
reactor  (PrOX)  for  the  gas  purification,  the  fuel  cell  which 
produces  electrical  energy  and  the  catalytic  burner  for  the 
conversion  of  residual  hydrogen  together  with  the  heat  ex¬ 
changers  necessary  for  thermal  coupling. 

With  respect  to  the  above  mentioned  requirements  there 
is  also  the  need  to  optimise  the  dynamic  behaviour  of  a  fuel 
cell  system  based  on  gasoline  reforming.  The  quicker  the 
system  is  able  to  follow  load  changes  the  smaller,  cheaper 
and  less  heavy  is  the  necessary  battery  in  the  car. 

For  that  purpose,  the  components  of  a  fuel  cell  system 
with  autothermal  reformer  have  to  be  described  by  appro¬ 
priate  dynamic  simulation  models.  These  models  have  to  be 
implemented  in  a  simulation  program  which  is  able  to  solve 
differential  equations  with  the  help  of  a  numeric  solver.  Af¬ 
ter  implementation  the  dynamic  reaction  of  the  modelled 
fuel  cell  system  on  load  changes  from  10  to  90%  of  full  load 
and  from  90  to  10%  of  full  load  are  calculated. 


2.  Development  of  dynamic  models  for  the  system 
components 

The  modelling  of  chemical  reactors  can  be  done  in  dif¬ 
ferent  ways  of  which  the  continuous  models  are  the  most 
important.  Continuous  models  can  be  classified  as  homoge¬ 
neous  and  heterogeneous  models.  Homogeneous  models  do 
not  distinguish  between  the  gas  phase  and  the  solid  (catalyst) 
phase  and  therefore  they  have  the  advantage  of  a  reduced 
number  of  numerical  equations  although  this  type  of  model 
is  not  suitable  to  locate  Hot  Spots  in  the  catalyst  phase  [9] . 
In  heterogeneous  models,  gas  phase  and  catalyst  phase  are 
modelled  separately  and  therefore  they  are  able  to  describe 
the  proceeding  of  a  chemical  reaction  in  detail  and  addition¬ 
ally  they  are  able  to  localise  and  quantify  Hot  Spots.  The  dis¬ 
advantage  is  the  high  number  of  numerical  equations  which 
may  lead  to  numerical  problems  during  the  simulation  [9]. 

Dynamic  models  for  the  system  components  are  devel¬ 
oped  in  the  following  manner: 

•  Because  the  reactions  which  occur  in  the  water  gas  shift 

reactors  are  not  very  complex  and  their  kinetics  are  known, 
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Fig.  1.  Basic  design  of  a  fuel  cell  system  with  gasoline  reformer. 


these  reactors  are  described  by  the  help  of  heterogeneous 
models. 

•  The  autothermal  reactor  is  modelled  by  the  help  of  a 
homogeneous  model  because  the  reactions  occurring  for 
which  simplified  kinetic  data  are  available,  are  very  com¬ 
plex. 

•  Considering  the  modelling  of  the  other  chemical  reactors 
no  reaction  kinetic  is  available  (PrOX,  Catalytic  Burner). 
In  case  of  the  fuel  cell,  the  process  description  would  re¬ 
quire  complex  models.  These  components  are  described 
by  the  help  of  ID  dynamic  models  based  on  integral  bal¬ 
ances. 

•  The  modelling  of  the  heat  exchangers  is  done  by  dynamic 
models  which  contain  partial  differential  equations  for  the 
hot  and  the  cold  medium  as  well  as  for  the  heat  transfer 
wall. 

In  the  following,  the  modelling  of  a  chemical  reactor  is 

described  as  an  example  for  the  water  gas  shift  reactors. 
The  assumptions  and  simplifications  made  are: 

•  The  pressure  in  the  reactor  is  considered  as  constant. 

•  The  gaseous  phase  is  composed  of  H2,  CO,  CO2,  N2,  H2O 
and  CH4  and  is  with  respect  to  the  low  pressure  and  the 
high  temperature  considered  as  ideal  gas. 

•  The  influence  of  radiation  is  neglected. 

•  The  properties  of  the  solid  phase  are  considered  as  con¬ 
stant. 

•  The  properties  of  the  gaseous  phase  depend  on  composi¬ 
tion,  pressure  and  temperature  and  are  described  by  the 
thermodynamic  properties  of  the  individual  components 
in  the  ideal-gas  state.  The  heat  conductivity  and  the  dy¬ 
namic  viscosity  are  approximated  by  the  properties  of  N2. 

•  The  reactor  is  considered  as  an  adiabatic  device. 

•  The  distribution  of  the  gas  on  all  channels  of  the  monolith 
is  ideal. 

•  There  are  no  secondary  reactions  implemented. 


Fig.  2  shows  a  balance  element  in  a  reactor  where  the 
catalyst  is  a  solid  phase  and  fixed  on  a  monolithic  carrier  with 
many  channels  parallel  to  each  other.  The  model  considered 
is  a  dynamic,  ID  model. 

The  following  equations  can  be  derived  from  the  balances 
for  energy  and  mass  which  are  similar  to  those  known  from 
other  authors  e.g.  [10]: 

Energy  balance  for  the  gaseous  phase: 

dTg  1  d  d2Tg 

epgc^g—  =  -j-(ngcp,gTg)+skg^ 

T  Ois,gHy{Ts  —  7g) 

Energy  balance  for  the  catalyst  phase: 


+(i-s)J2rj^hj 

j 


Mass  balances  for  the  gaseous  phase: 


£Pg 


fyi,  g 

dt 


1  9  . 

~  ~A~dx 

—  PiPgCLyiyi^g  —  )7,  s) 


Ty 


K,xiK,x 

A 

i 

4>y 

K,x+dxiK,x+dx 
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Fig.  2.  Reactor  element  for  the  shift  reactor  models. 
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Mass  balance  for  the  catalyst  phase: 

(1  —  s)Pg  ^  =  fiiPgdy  (y/,g  —  yi, s)  —  (1  —  s)  y  ]  ( vi,jRj ) 

j 

As  start  conditions  for  t  =  0  and  0  <  x  <  L  profiles  for  the 
variables  in  the  reactor  are  implemented: 

Tg(x,  0)  =  7gV),  Ts(x,  0)  =  Tfx). 

yi, g(x,  0)  =  ylg(x )  yu six,  0)  =  yfs(x) 

For  the  entrance  of  the  reactor  boundary  conditions  after 
Danckwerts  [11]  are  implemented. 

The  kinetics  for  the  water  gas  shift  reaction  are  derived 
from  experimental  results  and  implemented  in  the  reactor 
models. 

The  implementation  of  the  reactor  models  in  the  sim¬ 
ulation  program  Matlab  Simulink  is  done  a  so-called 
S-function.  The  differential  equations  are  solved  with  a 
numeric  algorithm  which  is  not  able  to  solve  partial  differ¬ 
ential  equations.  For  that  reason,  a  numerical  solution  of 
the  partial  differential  equations  is  required.  The  numerical 
solution  is  done  in  the  variable  of  space  by  the  help  of  the 
method  of  finite  differences. 

3.  Results  and  discussion 

In  the  following,  the  results  obtained  for  a  thermally  in¬ 
tegrated  system  which  is  shown  in  Fig.  3  are  presented. 
With  the  developed  system  models  which  are  implemented 
in  the  simulation  program  Matlab  Simulink  load  changes 
from  10  to  90%  of  full  load  and  from  90  to  10%  of  full 
load  are  calculated.  In  this  article,  as  an  example  the  re¬ 
sults  for  a  load  change  from  10  to  90%  of  full  load  will  be 
presented. 


Fig.  4.  Reformate  flow  at  a  load  change  from  10  to  90%  of  full  load. 


Fig.  5.  Educt  flow  at  a  load  change  from  10  to  90%  of  full  load. 


Fig.  4  shows  the  reformate  flow  at  a  load  change  from 
10  to  90%  of  full  load.  The  new  steady-state  after  the  fuel 
cell  is  reached  after  approximately  20  s.  The  reason  for  this 
behaviour  is  made  clear  by  Fig.  5  which  shows  the  educt 
flow  at  this  load  change.  The  educt  flow  at  the  outlet  of 


Air 


Fig.  3.  Thermally  integrated  system  used  in  the  calculations  of  load  changes. 
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Fig.  6.  Vapour  number  at  a  load  change  from  10  to  90%  of  full  load. 


Fig.  8.  Reformate  composition  at  ATR-outlet  at  a  load  change  from  10 
to  90%  of  full  load. 


the  heat  exchanger  HX-1  rises  for  the  first  time  when  the 
higher  air  flow  out  of  the  compressor  K- 1  has  flown  through 
HX-1.  The  response  time  of  the  educt  flow  is  dominated 
by  the  residence  time  for  the  demanded  water  in  the  heat 
exchangers,  which  is  mainly  influenced  by  heat  exchanger 
HX-5  which  is  a  large  heat  exchanger  and  heat  exchanger 
HX-2  where  the  water  is  nearly  not  vaporised.  So  also  the 
reformate  flow  is  influenced  by  the  educt  flow  and  shows  a 
similar  behaviour. 

Fig.  6  shows  the  vapour  number  of  the  educt  flow  (water) 
at  the  load  change.  The  vapour  number  is  defined  as  water 
enthalpy  standardised  on  the  state  of  boiling: 

h  —  h' 

vapour  number  =  - . 

Ahb 

The  educt  water  is  fully  vaporised  after  it  has  flown  through 
the  heat  exchanger  HX-5.  In  the  heat  exchanger  HX-2,  it  is 
nearly  liquid  and  this  causes  a  high  residence  time  because  of 
the  high  density.  The  vapour  number  reaches  its  steady- state 
after  approximately  50  s. 

Fig.  7  shows  the  educt  temperature  at  the  outlet  of  heat 
exchanger  HX- 1  and  the  inlet  of  the  ATR  and  the  reformate 
temperature  at  the  ATR  outlet.  The  temperature  of  educt  at 
the  outlet  of  the  heat  exchanger  HX-1  only  shows  a  slight 
and  soft  increase.  The  temperature  of  educt  at  the  ATR  inlet 


is  influenced  by  the  different  residence  times  for  air  and 
water  in  the  system.  At  a  load  change  from  10  to  90%  of 
full  load  gasoline  is  mixed  with  a  lower  stream  of  air  and 
water  which  causes  a  decrease  in  temperature.  It  reaches 
steady-state  after  25  s  when  the  flow  of  air  and  water  has 
reached  their  normal  level.  In  consequence,  the  operation 
parameters  for  the  ATR  are  temporarily  shifted  to  abnormal 
conditions  (very  low  S/C,  normal  SR).  This  is  the  reason 
for  the  increase  of  the  ATR  outlet  temperature  in  the  range 
of  100  K.  In  addition,  temperature  peaks  especially  at  the 
catalyst  front  where  partial  oxidation  of  gasoline  is  the  main 
reaction  are  possible. 

The  consequences  on  the  product  composition  at  the  ATR 
outlet  are  shown  in  Fig.  8.  The  mole  fractions  for  CH4  and 
CO  increase,  whereas  the  H2  mole  fraction  decreases.  The 
short  peak  for  CH4-mole  fraction  is  a  consequence  of  the 
low  SR  in  the  first  2  s  after  the  load  change.  Because  of 
the  low  S/C  in  the  first  20  s  after  the  load  change  the  water 
gas  shift  reaction  in  the  ATR  is  obstructed  and  consequently 
less  CO  is  converted  into  H2.  A  normal  product  composition 
is  reached  after  25  s  when  the  operating  parameters  have 
reached  the  demanded  conditions. 

Fig.  9  shows  the  temperature  and  the  CO  concentration 
in  the  shift  reactors  at  a  load  change  from  10  to  90%  of 


a 

o 


o 

o3 


I 


s 


-■-THTS  -*-T  LTS  ^>CO  HTS  ^CO  LTS 


Fig.  7.  Temperature  at  HX-1 -outlet,  ATR-inlet  and  ATR-outlet  at  a  load  Fig.  9.  Shift-reactor  temperature  and  CO-concentration  at  a  load  change 
change  from  10  to  90%  of  full  load.  from  10  to  90%  of  full  load. 
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Fig.  10.  Electrical  power  from  fuel  cell  at  a  load  change  from  10  to  90% 
and  from  90  to  10%  of  full  load. 

full  load.  The  temperatures  of  the  shift  reactors  are  mainly 
influenced  by  the  performance  and  the  thermal  mass  of 
the  heat  exchangers  placed  in  front  of  them  (HX-l/HX-2). 
For  that  reason,  they  reach  steady- state  after  approximately 
300  s.  Because  of  the  lack  of  water  in  the  heat  exchanger 
HX-1  during  the  first  20  s  the  reformate  cannot  be  cooled 
down  so  far  as  foreseen  in  the  design  state,  consequently  the 
temperature  of  the  HTS  rises.  In  the  heat  exchanger  HX-2, 
the  whole  amount  of  water  is  available  immediately,  the  re¬ 
formate  is  cooled  too  strongly.  The  temperature  in  the  LTS 
decreases.  Due  to  the  high  CO  production  in  the  ATR  also 
the  CO  concentration  at  the  outlet  of  the  shift  reactors  is 
much  higher  than  in  the  design  state.  Steady-state  for  the  CO 
concentrations  in  these  reactors  is  reached  after  nearly  100s. 

Fig.  10  shows  the  electrical  power  available  from  the  fuel 
cell  for  the  two  load  changes.  The  electrical  power  available 
from  the  fuel  cell  depends  on  the  reformate  mole-flow  and 
the  H2  concentration.  For  a  load  change  from  10  to  90%, 
it  reaches  steady-state  after  20  s,  for  a  load  change  from  90 
to  10%  it  is  reached  after  3  s.  The  curves  for  the  reformate 
mole  flow  are  reflected.  The  curve  for  the  load  change  from 
10  to  90%  also  reflects  the  decrease  of  H2  concentration  due 
to  the  too  low  SR  and  S/C  in  the  ATR. 

The  results  for  a  load  change  from  90  to  10%  of  full  load 
are  based  on  the  same  effects,  but  the  response  times  are 
longer  because  of  the  lower  gas  flows. 

4.  Conclusion 

According  to  the  results  obtained  the  dynamic  behaviour 
of  the  reformer  is  mainly  influenced  by  the  volume  of  the 
different  components  which  has  an  influence  on  the  resi¬ 
dence  times.  The  calculations  also  show  that  the  dynamic 
behaviour  is  also  influenced  by  the  mass  of  each  compo¬ 
nent  whose  high  heat  capacity  on  the  one  hand  leads  to  a 
function  as  an  energy  storage  and  on  the  other  hand  absorbs 
changes  in  temperature  with  the  result  that  it  takes  much 
longer  for  the  temperature  to  reach  the  new  steady- state  but 


also  prevents  the  component  from  temperature  peaks  as  it 
works  like  a  buffer. 

In  order  to  make  the  reformer  as  dynamic  as  possible  one 
has  to  construct  a  compact  and  light-weight  reactor  system 
which  enables  short  residence  times  of  the  gases  in  the  dif¬ 
ferent  sections  and  short  times  for  the  temperature  to  reach 
its  new  steady- state. 

The  calculated  response  times  also  show  that  even  in  re¬ 
gard  of  the  action  items  described  in  order  to  construct  a 
dynamic  reformer  the  dynamic  behaviour  of  the  system  will 
be  too  slow  to  work  as  a  propulsion  system  without  any  en¬ 
ergy  storage.  The  required  system  dynamic  for  mobile  ap¬ 
plications  is  a  response  time  of  1  s  for  a  load  change  from 
10  to  90%  of  full  load.  So  it  can  be  easily  estimated  that  the 
required  dynamic  of  such  a  system  will  only  be  achieved 
with  the  help  of  a  battery  which  is  electrically  loaded  by 
the  fuel  cell  system.  Nevertheless  it  is  absolutely  necessary 
to  build  a  fuel  cell  system  which  is  as  dynamic  as  possible 
because  with  increasing  dynamic  of  the  fuel  cell  system  the 
required  capacity  and  power  of  the  battery  becomes  lower 
and  the  battery  itself  becomes  cheaper  and  lighter. 
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